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ABSTRACT 

We present the prompt BAT and afterglow XRT data of 5w//f-discovered GRB 08 1 109 A up to 
~ 5 x 10 5 sec after the trigger, and the early ground-based optical follow-ups. The temporal 
and spectral indices of the X-ray afterglow emission change remarkably. We interpret this 
as the GRB jet first traversing the freely expanding supersonic stellar wind of the progenitor 
with density varying as p oc r~~ 2 . Then after approximately 300 seconds the jet traverses into 
a region of apparent constant density similar to that expected in the stalled-wind region of a 
stellar wind bubble or the interstellar medium (ISM). The optical afterglow data are generally 
consistent with such a scenario. Our best numerical model has a wind density parameter A* ~ 
0.02, a density of the stalled wind n ~ 0.12 cm -3 , and a transition radius ~ 4.5 x 10 17 cm. 
Such a transition radius is smaller than that predicted by numerical simulations of the stellar 
wind bubbles and may be due to a rapidly evolving wind of the progenitor close to the time 
of its core-collapse. 

Key words: Gamma Rays: bursts -GRBs: individual (GRB 081 109A)-ISM: jets and 
outflows-radiation mechanisms: nonthermal 



1 INTRODUCTION 

Gamma-ray bursts (GRBs) are the most luminous explosions in the 
Universe. They feature extremely relativistic outflows with bulk 
Lorentz factors of 10 2 ~ 3 and isotropic energies of 10 48 ~ 54 ergs. 
These energies are widely believed to be generated via the core- 
collap se of massive s tars for conventional long (-duration) GRBs 
(e.g.. IWooslev1ll993l : Iwooslev & Bloomll2006h or the merger of 
compact star binarie s for conventional short(-duratio n) GRBs (e.g., 
lEichler et all 1 19891 : iNaravan. Paczvnski. & Prianl 1 19921) . In the 
standard fireball model the prompt soft 7-ray emission is pow- 
ered by the collision of the material shells within the relativis- 
tic outflow (i.e., the internal shocks); afterwards these material 
shells spread and merge into a single uniform outflow, contin- 
uing to move outwards; then the long-lived X-ray, optical and 
radio afterglow emission is powered by the intera ction of the 
overall outflow and the circum-burst med ium (e.g.. |Piran1ll999l : 
Meszaros 2002; Zhang & Meszaros 2004). Therefore, the tempo- 
ral and spectral evolution of the multi-wavelength afterglow can be 
used to diagnose the underlying radiation mechanism and the pro- 
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file of the circumburst medium (e.g., Sari, Piran, & Naravan 1 19981: 
IChevalier & Lill2000l:|Panaitescu & Kumar l2002h . 

For short GRBs, the circum-burst medium is expected to be 
of the interstellar medium (ISM) type or even of the intergalac- 
tic medium (IGM) type, for which the number density is roughly 
constant and much lower than unity, being consistent with the 
current afterglow modeling (see Nakarll2007l for a review). For 
long GRBs, Panaitescu & Kumar (2002) found that the after- 
glow modeling usually favors the constant density (CD) medium 
scenario, half cases in 10 bursts of a sample were better fitted 
by the CD medium while only one better described by a free 
wind (FW) medium. However, Starling et al. (2008) found that 
in 10 bursts of a sample, 4 were clearly in FW medium while 
only 1 consistent with CD medium. It seems that the circum- 
burst environments are unlikely to be drawn from only one of the 
CD or FW profile for long bursts. This is not obviously consis- 
tent with the collapsar model in which an ideal FW circum-burst 
medium would be created. The FW has a density profile p — 
5 x Kr^r-'gcm-'.teA = ( TJJ ^_ T )( 1222^1) 
is the parameter reflecting the density of the wind. One poten- 
tial solution is that the ideal FW medium profile has been mod- 
ified before the GRB explosion. As is known, massive stars are 
believed to enter the Wolf-Rayet stage during their late evolu- 
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tion and have lost a major fraction of their masses in the form of 

the stellar wind. The interaction between this stellar wind and the 

I 1 

surrounding medium creates a bubble structure (e.g. Weaver et al. 



Chevalier, Li & Fransson 




2004 Ivan Marie et al. 1 l200d 2001 


2008;Eldridgeetal. 2006; 


Eldridge 2007; Pe'er & Wiiers 12006). 



In this scenario the free expanding supersonic wind is terminated 
at a radius R t ~ 10 18 - 10 2() cm, where the density jumps by a 
factor of 4 or more with the lower value expected for such an adi- 
abatic shock. Beyond the wind termination shock is the roughly 
constant density stalled wind, holding up to a rather large radius 
Rism ~ 10 19 -10 21 cm, outside this radius is the very dense swept- 
up ISM (the number density n ~ 10 2 — 10 3 cm -3 ) and then the 
ISM. 

If the above picture is a good approximation of the circum- 
burst medium, it may be possible to observed afterglow signatures 
caused by the medium transition at ~ R t or even at ~ -Rism- In 
Section 2, we describe the signatures in GRB afterglows for the 
first transition which is more easily observed and more of interest 
with respect to the second one. Such a transition is possibly evident 
in GRB 050319 (Kamble et al. 2007) and GRB 050904 (Gendre et 
al. 2007). In Section 3, we present the X-ray afterglow observations 
of GRB 081 109 A which leads to reach the conclusion that the FW- 
CD transition is quite clear in this event. The optical afterglow data 
are generally consistent with such a scenario. Throughout this work 
we use the notation F(t, v) oc t a v® for the afterglow monochro- 
matic flux as a function of time, where v represents the observer's 
frequency, a is the monochromatic flux decay index, and /3 is the 
energy spectrum index. The convention Q x — Q/10 x has been 
adopted in cgs units. 
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Table 1. Temporal index a and spectral index /3 in the cases that the 
medium profile can be constrained, where the convention F(t, v) oc 
t a v@ is adopted, p is the power-law index of the swept-up electrons, and 
v a , v c , and v m are the self-absorption frequency, the cooling frequency 
and the characteristic frequency in the synchrotron radiation, respectively 
iSari, Piran, & Naravanll 998; Chev alier & Li |2000|) . 



Case 

V & < V < V c < V m 
V a < V < V m < V c 
max{^ a , Vm} < V < v c 

Vc < V 



tion at the radius Rt ? It is known that the characteristic frequency 
Vm decays with time as i~ 3//2 in both CD and FW scenarios. How- 
ever, the cooling frequency v c evolves very differently. In the FW 
scenario it increases as t 1 ^ 2 but in the CD scenario it declines as 
£-1/2 Therefore, the following observational signatures would be 
evident before and after the transition. In Case I, the deceline of 
the X-ray emission remains unchanged while the optical decay be- 
comes shallower by a factor of t 1 ^ 2 as long as the optical band is 
still below v c . In Case II, the X-ray decay will get flattened by a fac- 
tor of t 1 ' 2 as long as the X-ray band is still below v c . Afterwards, 
when v c drops below the X-ray band, the X-ray decay steepens by a 
factor of t~ 1/l4 . As we'll show in next Section, the X-ray afterglow 
of GRB 081 109 A fits Case II. 



2 THE AFTERGLOW SIGNATURES OF THE MEDIUM 
TRANSITION 

The GRB afterglow emission in the CD or FW scenarios has been 
extensively discussed (see Zhang & Meszaros 2004 for a review). 
We use the standard fireball afterglow theory with the simple mi- 
crophysical assumptions of constant energy fractions imparted to 
the swept-up electrons, e c , and to the generated magnetic field, 
£b, respectively. The typical synchrotron radiation frequency v m 
and the cooling frequency v c are calculated in the standard way 
dSari, Piran. & Naravan|[l998l) . Both scenarios lead to the after- 
glow closure relations made of the temporal decay index a and the 
spectral index j3, depending upon the spectral segment and the en- 
ergy distribution index of p ~ 2 — 3 of the swept-up electrons. Usu- 
ally the synchrotron self-absorption frequency, V&, is much lower 
than the optical band and thus neglected unless radio observation 
invoked. Therefore, with good-quality temporal and spectral data, 
the circum-burst medium profile can be reliably constrained. We 
summarize the closure relation in Ta ble [T] apart from the v < v a 
cases (see also lPe'er & Wiiers1l2006t) . 

For the X-ray afterglow of our interest, there are two ways to 
diagnose the number density profile of the medium: 

(I) When v 1T1 < v op t < v c < fx occurs, in the FW scenario 
the X-ray emission drops with time as t( 2_3p )/ 4 while the optical 
emission drops faster by a factor of 1/4. On the contrary, in the CD 
scenario the optical emission drops more slowly than the X-rays. 

(II) In the case of v m < vx < v c , i.e., the X-ray emission is in 
the slow cooling phase, the temporal and spectral indices roughly 
satisfy a — (3/3 — l)/2 in the FW scenario. 

Then how does the afterglow evolve after the FW-CD transi- 



3 THE FW-CD TRANSITION IN GRB 081109A 

3.1 Observation and data reduction of GRB 081109A 

GRB 081 109 A was triggered and located by Burst Alert Telescope 
(BAT) onboard S wift at 07:02:06 UT (trigger 334112), Novem- 
ber 9th, 2008 (see llmmler et al. 1120081) . We downloaded the raw 
data from the UK Swift data archive and processed them in a 
standard way with HEAsoft 6.6. The BAT lightcurve in the 15- 
350 keV band was processed with the batgrbproduct task. 
This burst has a Tgo about 61 seconds in 15-350 keV, classified 
as a long GRB. The duration value is consistent with the dura- 
tion measurement of 45 s in 8-1 000 keV by the Gamma -ray Burst 
Monitor (GBM) onboard Fermi jvon Kienlin et al. 120081) . Also the 
time-averaged GBM spectrum is best fit by a power-law function 
with an exponential high energy cutoff. The power law index is - 
1.28±0.09 and the cutoff energy, parameterized as E p is 240±60 
keV (chi squared 510 for 478 d.o.f.). The fluence in 8-1000keV is 
6.35 ±0.43 x 10 -6 erg-cm -2 . Fig.Q]shows the prompt lightcurves 
of GRB081109A in 15-25, 25-50, 50-100 keV with 1 s time bin- 
ning. Spectral lag between different energy bands cannot be well 
measured for this event due to relatively low signal to noise. 

The Swift X-ray Telescope (XRT) observations began at 65.6 
seconds after the BAT trigger and discovered a bright and fad- 
ing X-ray afterglow. Observations continued during the following 
hours and days and in several return visits, with Windowed Tim- 
ing (WT) mode for ~ 300 s after the trigger and Photon Counting 
(PC) mode afterwards. Throughout the X-ray observation, spec- 
tral softening is evident. There is no need for pile-up correction 
to WT data as the highest count rate is less than ~150 count s , 
but such correction should be applied to early PC data when the 
count rate is higher than ~ 0.6 count s _1 in order to get correct 



3 




10 20 30 40 50 

Time since trigger (sec) 



Figure 1. The prompt lightcurves of GRB 081 109A in 15-25, 25-50, 50-100 
keV with 1 s binning. Spectral lag between different energy bands cannot 
be well measured for this event due to relatively low signal to noise. 



X-ray lightcurve and spectra. We made this correction by fitting 
a King function profile to the point spread function (PSF) to de- 
termine the radial point at which the measured PSF deviates from 
the model. The counts were extracted using an annular aperture 
that excluded the affected ~4 pixel core of the PSF, and the count 
rate was corrected according to the model. We also considered the 
most recent calibration and exposure maps. Fig.[2] shows the 0.3- 
10 keV X-ray lightcurve, which can be well modeled with a dou- 
bly broken power-law which is equivalent with a broken power- 
law with a jump transition at ~500 s. For a doubly broken power- 
law, the fitted parameters are: ai = -1.75 ± 0.04 (x 2 = 152.4 
for 119), t b i ~ 310 s, a 2 = -0.70 ± 0.13 {\ 2 = 13.7 for 
18), i b2 ~ 2.9 x 10 3 s, a 3 = -1.24 ± 0.03 ( X 2 = 38.9 for 
54). The integrated spectra of the above first and third segments 
are shown in Fig.[3] In detail, the spectral power-law indices are 
Pi = -0.74 ± 0.05, p 3 = -1.27 ± 0.10, and /3 2 is in the middle. 

The Swift Ultraviolet and Optical Telescope (UVOT) began 
observing at ~150 s after the trigger and fo und no optical coun- 
terpart down to ~18 mag jlmmler et al. 1120081) . Ground-based ob- 
servations of the afterglow of GRB 081 109 A were carried out with 
the REM telescope at La Silla l lZerbi et ail 1200 ll ; Ichincarini et all 
120031 ICovino et al. ||2"004» equipped with the ROSS optical spec- 
trograph/imager and the REMIR near-infrared camera on 2008 
Nov 9, starting about 52 seconds after the burst jD' Avanzo et al. I 
2008). The night was clear, with a seeing of about 2.0". We col- 
lected images with typical exposures times from 10 to 120 sec- 
onds, covering a time interval of about 0.5 hours. Image reduction 
was carried out by following the standard procedures. Astrome- 
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Figure 2. The X-ray afterglow lightcurve of GRB 081 109A in 0.3-10 keV. 
Also marked are the fitting temporal indices for the WT and late PC seg- 
ments of the lightcurve. 
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Figure 3. The spectra of the WT and late PC segments of the X-ray af- 
terglow of GRB081109A. An absorbed power-law fitting is adopted. Also 
marked are the fitting spectral indices respectively. 
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Table 2. Optical photometry for GRB 081109A. Upper limits are at 3a 
confidence level. 



Filter Time since trigger (s) Time integral (s) Magnitude 



H 


89.4 


71.3 


>15.0 


H 


168.9 


72.2 


15.017±0.115 


H 


248.4 


71.3 


>15.4 
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327.0 


71.3 


14.893±0.109 
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394.4 


47.3 


14.747 ±0.128 
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71.3 


14.963±0.114 


H 


990.6 


121.1 


>15.5 


H 


1545.3 


171.7 


15.763±0.101 


K 


450.6 


45.6 


14.365±0.259 


K 


718.4 


72.2 


14.465±0.207 


K 


1122.8 


121.1 


14.645±0.188 



try was performed using the USNOBl.cE and the 2MAS50 cata- 
logues. We performed aperture photometry for the afterglow and 
comparison stars. The afterglow was not detected in the optical. In 
the NIR it was detected only in the H and Ka bands. Results for 
the photometry are reported in Table|2] The GROND instrument 
equipping the 2.2m ESO/MPI telescope at La Silla started obser- 
vation 17.1 hr after the trigger, gave a redshift limit of z < 3.5 
and a best fit of intrinsic extinction of Av between 0.6 and 1.2 
(CI emens. Kruehler. & Greiner 2008). 



3.2 Analytical investigation of the afterglow of GRB 081109A 

For t < tbi, the temporal index on = —1.75 ± 0.04 and the 
spectral index (3i = -0.74 ± 0.05 are related by a - 3/3/2 w 
—0.5, being consistent with the forward-shock emission in the FW 
medium as long as i/ m < f x < v c . For t > tb2, the temporal index 
q 3 = -1.24 ± 0.03 and the spectral index /3 3 = -1.27 ± 0.10 
are related by a — 3/3/2 ~ 0.5, suggesting tha t the medium can be 
either FW or CD given u x > max{f m , is c } (Zhang & Meszaros 
l2004h . If the FW scenario holds from the very beginning, we have 
v c {th2) > t'c(^bi) > ^ x because of v c oc i 1 / 2 . One can, of course, 
assume that either y4* or e B has increased abruptly and then get a 
J / C (^b2) <C v c {tbi)- However, such a treatment is lack of any solid 
physical background and is thus artificial. On the other hand, if we 
assume that at t ~ tbi there comes the FW to CD transition, the 
forward shock emission light curve will get flattened by a factor 
of t 1//2 , roughly consistent with the data (please note that Q2 is 
only poorly constrained). In this scenario, the steepening at t ^ 
th2 implies that v c (tb2) > v x - Below we present our quantitative 
estimates. 

In a FW medium, we have ai = p and (3i — — in the 
case of v m < ux < Vq- In a CD medium, we have «3 = and 
/?3 = — | for i/x > max{^ m , v c }. We find that p = 2.5 fits both 
the temporal and spectral slopes of GRB 81 109A. 

In a FW medium, we have (e.g.. lChevalier & Li 12 000): 

-Fi/.max = 0.23 Jy D L 2S 34 e-J_ 2 Ej s3 A t t d / 3 , (1) 

, Q „ in 16 Ti ^,2^1 + ^,1/2 2 1/2 ,-,1/2.-3/2 

v m = 1.8x10 HzC p (— — ) e e ,-ie B ,-2^,53^,-3 > ( 2 ) 



1 http://www.nofs.navy.mil/data/fchpix/ 

2 http://www.ipac.caltech.edu/2mass/ 



l 1 „ ml3 u t 1 "T" Z \-3/2 -3/2 j-,1/2 ,-2,1/2 , Q1 
U C = 1.1 X 10 Hz (— — ) e B ,-2- B k,53 A » *d'-3> ( 3 ) 

where is the isotropic equivalent energy, z is the redshift of 
the GRB and _Dl is the corresponding luminosity distance, td is 
the time in days since trigger in the observer's frame, and C v = 
13(p - 2)/[3(p - 1)]. 

In our model, at t ~ 65.6 s, i/ m < 0.3keV, v c > lOkeV and 
fb.3keV ~ 2.6 x 10" 3 Jy are needed. So we have 

eo,-lEB,-2£k,53 < 3 (Um < 0.3keV), (4) 

t^E'l^A* < 0.002 (i/ > lOkeV), (5) 

^-i^-2<53^* ~ 0.01. (6) 

In a CD medium, we have (e.g., ISari. Piran. & Naravan I 

Il998h : 

-F^max = 6.6 mjy (^-^-)D^ 2 2S . S4 .^^ 2 E kt5 3nl /2 , (7) 

„ / „ 1n l6 n ^2/1 + 2x1/2 2 1/2 „l/2 .-3/2 
V m = 2.4X10 Hz C p (— — ) £ c ,-l£ B -2-^,53^, -3 > ( g ) 

, , „ 1n l6u /l + Z\-l/2 -3/2 d-1/2 -1.-1/2 /n , 

u c = 4.4x10 Hz(— ^— ) e B ,-2 B k,53 n o *d -3. ( 9 ) 

where n is the density of CD medium. 

At t ~ tb2 ~ 2900s, our model suggests that i/ m and < 

0. 3keV, and F .3 kcV ~ 5.7 x 10" 5 Jy. We then have 

4,-i£B,- 2 £k,53 < 8 x 10 4 (v m < O.SkeV), (10) 
e B ,-2#k,53rJo > 0.01 (^ c <0.3keV), (11) 

3/2 1/8 p 9/8 , a 
e o _i£ B -2^,53 ~ L -°- ( - lI > 

In a termination shock model, the crossing time is estimated 
bv lChevalier. Li & Fransson I d2004h : 

t(R t ) = 1.5h(ii£) B^A^mo 1 ~ 310s, (13) 

1. e., 

S k J 3 A 2 ,_mo 1 ~ 0.06. (14) 

For GRB 08 1109 A, there is no self-consistent solution for 
Eqs. J4l6b . Eqs. d 10112b and Eq.J14b provided that e B is a constant in 
the free wind and in the CD medium. With eqs. CD, l|5) and fill, 
we have £b,cd /e B ,w ^ 7, where the subscripts "CD" and "w" rep- 
resent the physical parameters measured in CD and FW medium, 
respectively. A similar assumption was needed in the mod eling 
of the aftergl ow data of GRB 05 0904 dGendre et al71 12007|) and 
GRB050319 dKamble et al~ll2007l) . The physical reason is that the 
CD medium has been heated by the termination reverse shock and 
then may be weakly magnetized. 

The optical data, though rare comparing with the X-ray ones, 
provide us a reliable test of the current afterglow model. As shown 
in Fig.[4] the H-band lightcurve is distinguished by a flat at early 
time and then a re-brightening at t ~ 400 s. These features are 
generally consistent with our FW to CD medium model. In the FW 
medium, a flat segment is expected if the observer's frequency v b s 
is above the synchrotron self-absorption frequency v & but below 
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(< Vc)- In the CD medium, a re-brightenin g is present if the re- 
lation i/ a < v i, s < u m < v c still holds (e.g. JZhang & Meszaros I 
2004). The re-brightening peaks when u m crosses the observer's 
frequency. After then the flux will drop with time as £ 3 ( 1_p >/ 4 , a lit- 
tle bit shallower than the simultaneous X-ray decline, as suggested 
by the H-band data. 



3.3 Numerical fit of the afterglow of GRB 081109A 

Th e code used here to fit th e X-ray light curves has been developed 
bv lYan.Wei&Fanl {2007), with small changes to adapt a density 
transition in surrounding medium. Assuming a FW to CD medium 
transition, we find out that the observation data can be reasonably 
reproduced with the following parameters (see Fig.|4]and Fig.[5]l: 
E k = 4 x 10 54 erg, the initial Lorentz factor 70 = 500, A, = 0.02, 
n = 0.12 cm" 3 , R t = 4.5 x 10 17 cm, e c = 0.02, p = 2.5, 
fB, w = 0.0002 and £b,cd = 0.001. The source is assumed to 
be at a redshift z = 1. As shown in Fig.[5] a sudden increase of 
ee for t ^ ibi is required to account for the jump of v c inferred 
from the X-ray data. A» ~ 0.01 and R t ~ 10 17 - 10 18 cm is 
lower than the typical value found in numerical simulations and 
may be due to a rapidly evo lving wind of the progenitor close to 
the time of its core- collapse (Eldridg eeTal. ll2006l : Eidridge 12001 
van Marie et al. ]Hooi). Assuming a redshift 8 = 1, the isotropic 
energy _E 7 in the energy range of 1 — 10000 keV is about 5 x 
10 52 ergs. The corresponding GRB efficiency E 1 /{E 1 + Ek) is 
~ 1%. Such a low efficiency, thoug h not typical, is still reasonable, 
as fou nd in previous estimates (e.g jFan & Piranll200d : ljm~& Fan 
120071) . 

The observed optical flux is lower than the value extrapolated 
from X-ray observation, which implies that there is extinction in 
optical band. To agree with the optical observation, our numerical 
fit for the H and Ks band has been adjusted by a factor of 0.3, 
which requires about 1 mag extinction caused by the GRB host 
galaxy at the observed H and Ks band. They are about 1 micron 
at the GRB host galaxy. It essentially means Eb-v ~ 0.8 — 0.9 
assuming Milky Way or SMC extinction curves. With a Galactic 
dust to gas ratio it would correspond to Nh ~ 5 x 10 21 cm -2 . The 
fit is shown in Fig.|4] The small differences around the termination 
radius may due to the sharp density jump we have assumed. 

Before the blast wav e reaches Rt, its radius can be estimated 
as jChevalier & Lill200ti) : 

R = 3.5 x l0 17 (i±i)- 1 / 2 J B^A7 1/2 ^ /2 cm. (15) 

For t(R t ) ~ 310 s, with eq.© we have R = R t ^ 10 22 [(1 + 
z)/2] ~ x ^ 2 A Z J 2 £^ 2 _ 2 cm. This implied that Rt depends on the un- 
determined redshift z weakly. A larger A» requires a smaller eb,w 
otherwise R t will be much smaller than a few x 10 17 cm, the lowest 
value expected in the numerical simulation. Since the derived eb,w 
is already as low as ~ 10~ 4 , a smaller value is less likely. That is 
why we will not consider the case of a FW parameter A» 2> 0.01. 

When a GRB jet finally enters into the very dense swept-ISM 
region, strong reverse sh ock may be formed and an afterglow re - 
brightening is expected jDai & Wu ||2003| ; IPe'er & Wiiers Il200d) . 
The X-ray observation for GRB081109A lasted to about 4.7 x 10 5 
s after the trigger and did not find obvious flux enhancement. In 
our numerical fit, the jet front reached a radius ~ 1.9 x 10 18 cm 
at such a late ti me, much smaller than Rism estimate d by numeri- 
cal simulations teldridge et al. l2O06fclmdridge l2007h . So the GRB 
outflow was still in the shocked wind material ejected in the Wolf- 
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Figure 5. Numerical fit to the spectral evolution in GRB081 109A afterglow. 
In the X-ray band(0.3-10keV), the spectrum between 65.6s and 310s is a 
single power law j3\ = 0.75, and after ~ 4800s is another single power 
law/3 3 = 1.25. 

Raye t stage with an approximately constant density j Weaver et al. I 

E973). 



4 CONCLUSION AND DISCUSSION 

Although some GRB s are inferred to occ u r in free stellar win d 
medium (e.g.. IPanaitescu & Kumar! I2002L IStarling et al~l |2008|) . 
most are inferred to occur in a medium with a constant num- 
ber density (even for some bursts associated with bright super- 
novae, see iFan 112008 . and the references therein), which may in- 
dicate that GRB outflows expand into the wind bubble rather than 
the ideal free stellar wind. As shown in Section 2 of this work, 
in some cases the X-ray afterglow data could shed light on the 
wind bubble structure. Therefore the X-ray afterglow observation 
since the early time is required to trace the profile of the circum- 
burst medium. However, in many Swift GRB events the early X- 
ray aft erglow deviates from the standard afterglow model signifi- 
cantly dNouseket ai~l f2006; Zha ng et al. 1120061) . For instance, the 
prolonged activity of GRB central engines would generate ener- 
getic X-ray fl ares that have outshone the regular forward shock 
emiss ion (e.g jFan & Weil2005l : lzhang et al. l2006r . lNousek et al. I 
2006). Fortunately, in GRB 081 109 A there is no flare accompany- 
ing the early X-ray afterglow. The temporal and spectral evolutions 
of the X-ray afterglow imply a medium transition at the radius 
Rt ~ 4.5 x 10 17 cm (see Section 3). Such a small R t implies 
a small wind parameter A* and a large p/k since 

Rt = 5.7 x 10 17 ( J w , )( h P/ \ 7 7 )- 1/2 Al /2 2 cm, 
v 10 3 kms- lA 10 6 cm 3 K ; *~ 2 

where p is the pressure in the shocked win d and k is the Boltzmann 

constant (Cheva lier. Li & Franssonll2004l) . Indeed A* is found to 

be as small as 10~ in our numerical fit. GRB081109A is thus a 

good candidate of long GRBs born in wind bubble. 

The rising behavior of the very early afterglow lightcurves can 

play an important role in probing the density profile of the medium. 

This is particularly the c ase if the reverse s hock optical emission i s 

very weak. As shown in I Jin & Fan! ( 120071) and lXue et al~1 ^2009), 

in the free wind scenario the very early optical rise is usually 

not expected to be faster than t 1 ' 2 while in the constant density 
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Figure 4. Numerical fit to the afterglow of GRB081 109A. The solid circles and squares are X-ray observations, and the solid line is our numerical fit to X-ray. 
The circles on the top are H band (empty) and Ks band (solid) observations, the triangles are upper limits in H band. The dashed and dotted lines are numerical 
fit to H band and Ks band data, respectively. 



medium the rise can be faster than t . Therefore the ~ t -like 
rise in the early optical/infrared/X-ra y afterglows of GRB 060418, 
GRB 060607A dMolinari et al. |2007|), GRB 060801 GRB 060926, 
GRB 0803 19C, and GRB 0804 13B dXue et al. 1120091) rules out the 
free wind medium for R > 10 16 cm. Thus the absence of the free 
wind signature in some GRB afterglows is still a puzzle. One possi- 
ble solution is that the mass loss rate due to the stellar wind before 
massive stars collapse is not a constant and might be much lo wer 
than previously assumed jEldridge et al. ll2006l : lEldridge1l2007l) . 
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